Non-ideal magnetohydrodynamic simulations of subcritical pre-stellar cores with non-equilibrium chemistry

Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and forming stars. Here, we modify publicly available version adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey 2008) include a detailed treatment non-ideal MHD study such collapsing prestellar cores. We implement two very extended non-equilibrium chemical networks, largest which is comprised $\sim$ 300 species includes description deuterium chemistry. The ambipolar-diffusion, Ohmic Hall resistivities then self-consistently calculated from abundances charged species. present series 2-dimensional axisymmetric simulations where vary model, cosmic-ray ionization rate, grain distribution. benchmark our implementation against ideal previously-published results. show that, at high densities ($n_{\rm{H_2}}>~10^6~\rm{cm^{-3}}$), ion that carries most perpendicular parallel conductivities not $\rm{H_3^+}$ as was previously thought, but instead $\rm{D_3^+}$.